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Abstract. In the first part of this paper, the different distinguishable pathways and 
regions of the single and sequential double ionization are determined and discussed. It 
is shown that there are two distinguishable pathways for the single ionization and four 
distinct pathways for the sequential double ionization. It is also shown that there are 
two and three different regions of space which are related to the single and double ion- 
ization respectively. In the second part of the paper, the time dependent Schrodinger 
and Newton equations are solved simultaneously for the electrons and the nuclei of 
H2 respectively. The electrons and nuclei dynamics are separated on the base of the 
adiabatic approximation. The soft-core potential is used to model the electrostatic 
interaction between the electrons and the nuclei. A variety of wavelengths (390 nm, 
532 nm and 780 nm) and intensities ( 5 x 10 14 Wcmr 2 and 5 x 10 15 Wcm~ 2 ) of the 
ultrashort intense laser pulses with a sinus second order envelope function are used. 
The behaviour of the time dependent classical nuclear dynamics in the absence and 
present of the laser field are investigated and compared. In the absence of the laser 
eld, there are three distinct sections for the nuclear dynamics on the electronic ground 
state energy curve. The bond hardening phenomenon does not appear in this classical 
nuclear dynamics simulation. 
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I. Introduction 

Atoms and molecules exposed to intense laser pulses reveal a vast wealth of fascinating 
phenomena for example: single and double ionization [TJ, above-threshold ionization [2], 
charge resonance enhanced ionization [3J, dissociative-ionization [HE], above threshold 
dissociation |6J, bond softening and hardening [7J El El [TO] and high order harmonic 
generation [llj. Scientific research regarding these phenomena resulted in a revolution 
in ultrashort laser pulses and molecular science with many broad applications such as 
the control of the molecular processes [12], generation of a few cycles of femtosecond 
and attosecond pulses, the emergence of the attophysics [T3J, and time- resolved imaging 
of molecular dynamics and reactions [14J. 

Our knowledge about the interaction of atoms and molecules with ultra-short 
intense laser pulse has proceeded like essentially all other fields of atoms and molecules 
from simple systems such as Hydrogen atom and Hydrogen molecule to complex ones 
such as proteins. The first step in the chain of studying many-electron systems is 
investigation of the simplest two-electron systems, i.e. helium atom and H2 diatomic 
molecule [151 [TBI H7]- The response of many electron atoms or molecules to the 
pulse of the laser field is often described by assuming that only one electron is 
active and responsible for the emission; this is called the single active electron (SAE) 
approximation. In this approximation, other electrons are assumed to contribute in 
the dynamics through a static screening potential. However, for full description and 
understanding the behaviour of the two-electron systems, it is necessary to consider both 
electron simultaneously without SAE approximation. This work plays a central role in 
developing our understanding of the interaction of many-electron atoms and molecules 
with ultrashort intense laser fields [HI [16], [TTJ . For this purpose, we need to solve the 
time-dependent Schrodinger equation (TDSE) for the two-electron systems. This work 
is impossible with available current computing power. Therefore, at the present time 
unavoidably in many researches, modelling of the interaction of two-electron systems 
with the laser pulse is accomplished by the so-called soft-core Coulomb potentials which 
makes the numerical solution of TDSE possible [T8l IT9] [201 |2T1 122] . 

In this work, we consider the nuclei as classical particles. This model reduces the 
complexity of the problem and helps to show the details of the dynamics of the electrons 
without the necessity that we get involved in complexities of the dissociative-ionization 
process [U 1231 I2U I2S1 EH]- In this research, the indistinguishability concept and the 
symmetry properties between two electrons will be demonstrated with some details. 
The main focus of attention in this paper is about the details of the classical nuclear 
dynamics and in a separate article the details of the electrons dynamics are represented 

mi- 

The paper is organized as follows: In the Sec. II and III, the details of the numerical 
solution of the TDSE and simulation box are described. In the Sec. IV, the results of 
the simulations for the field-free case and also for different intensities and wavelengths 
are presented and discussed. Finally, the conclusion appears in Sec. IV. We use the 
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atomic units (h = m e = e = 1) throughout this paper unless it is stated. 
II. Numerical solution of the TDSE 

In this section we introduce the details of the numerical solution and also the simulation 
box that we used for investigating the dynamics of the hydrogen molecule in the presence 
of a strong laser field. In a linearly polarized laser pulse with intensity upto 1 x 10 17 
Wcm~ 2 , most of the electron and nuclei dynamics occur in the direction of the laser 
field [21]. Therefore, we choose a one-dimensional model for both the electrons and 
nuclei coordinates. In what follows, Ri and R2 indicate the nuclei positions and z\ 
and z 2 indicate the electron positions. Furthermore, M and m indicate the nuclei and 
electrons masses, respectively, and e is the electron charge. The temporal evolution of 
the electronic parts of such a system is described by the time dependent Schrodinger 
equation (TDSE) on the base adiabatic approximation, i.e. j2TJ 128] 

. di/)(z 1 ,z 2 ,t;R 1 (t),R 2 (t)) 
dt 

H e ( Zl ,z 2 , t; R 1 (t),R 2 (t))ij(z 1 ,z 2 , t; R^t), R 2 {t)) (1) 
where the electronic Hamiltonian for this system, H e (zx, z 2 , t; Ri(t), R 2 (t)), is given by 

H e (z l7 z 2l t; R l (t),R 2 (t)) = 

' d 2 d 2 ~ 



2m 



dz\ dz 2 



+ Vc(zx, Z2 ,t;R l (t),R 2 (t)), (2) 



V c (z u z 2 ,t; R^Riit)) 
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i,a=l 



Zi-R a ) 2 + aJ y(z!-z 2 ) 2 + b 
Z Z 

1 2 = + Oi + z 2 )E f(t) cos(wt), 
(Ri - R 2 ) 2 + c 



(3) 



where Z\ and Z 2 = 1 are the nuclei charge and a, b, and c are the soft-core parameters 
which are related to electron-nuclei, electron-electron, and nucleus-nucleus interaction, 
respectively . As this equatin shows, we use a soft-core potential to model the 
electrostatic interaction among the electrons and the nuclei by the screening parameters 
a, b and c. We choose a = b = 1.0, c = 0.03 [21]. The ground state energy and 
the equilibrium internuclear distance of H 2 become -1.39 and 2.13 respectively for this 
choice of the parameters (see Fig. [TJ. 

The laser-molecule interaction is formulated in the dipole approximation in length 
gauge (— er.E(t)), where E is the laser peak amplitude and oj = litv is its angular 
frequency. The envelope of the laser pulse, f(t), is set as 

f(t) = W(^tt) (4) 
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Figure 1. Calculated the electronic ground states surface of hydrogen molecule on 
the base a one-dimensional model for both the electrons. 

where T\ is the time duration of the field irradiation, set at ri=8 cycles in this work. 
1 cycle is equal to 1.30 fs ( 53.78 a.u.) and 1.77 fs ( 73.36 a.u.) for 390 nm and 
532 nm respectively. After the mentioned time, the simulation continues for 8 more 
cycles in which there is no laser field as shown in Fig. [2] to follow the behavior of the 
system after turning off the laser field. In this simulation, the time step is set to 5t= 
0.02. The differential operators in Eq. (2) are discretized by the 11-point difference 
formulas which have tenth-order accuracies [2H] • To solve the above TDSE numerically, 
we adopted a general nonlinear coordinate transformation for electronic coordinates. 
For the spatial discretization, we have constructed a finite difference scheme with a 
nonuniform (adaptive) grid for z\ and z 2 electronic coordinates, which are finest near 
the nuclei and coarsest at the border regions of the simulation box. More details of our 
calculations are described in our previous reports [2SJ ED]- The absorber regions are 
introduced by using fourth-order optical potentials at the zi and z 2 boundaries, in order 
to capture the photoelectrons and prevent the reflection of the outgoing wave packets 
at the borders of the grid. More details of our calculations are described in our previous 
work [30] . 

In this work, we separate the electrons and nuclei dynamics on the base of the 
adiabatic approximation [211 128] . We consider the electrons dynamics in the quantum 
approach while the nuclei dynamics are investigated in the classical manner. So, the 
time dependent Schrodinger and Newton equations are solved simultaneously for the 
electrons and the nuclei respectively. The equations govern the nuclear dynamics are as 
follows |2Tj: 

MR^t) = F 12 (t) + F le (t) + F L (t) 

MR 2 (t) = F 21 (t) + F 2e (t) + F L (t) (5) 
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Figure 2. The used laser electric field in this article has a sinus second order 
envelop function with 8 optical cycles. 



where F L (t) = eE f(t) cos(wt) is the laser forces exerted on each nucleus and F 2 i(t) 
—F\2(t) are the internuclear repulsion 

R 2 (t)-R 1 (t) 
(c+ | R 2 (t)-R 1 (t) 

and F ne is the attractive forces between the electrons and nuclei where e and n are 
related to the electrons and nuclei respectively 



F 21 (t) = - Fl2 (t) = - , ^{ m s/ 2 (6) 



FUt) = -JJ [a +[ R n (t)- Zl? r 12 
[R n (t)- Z2 ]\*( Zl ,z 2 ,t)\Z dzid ^ 



[a + [R n (t) - z 2 Yf 2 
These Newtonian equations are solved by the Verlet algorithm [3l] . 



III. Different pathways and regions of the single and sequential double 
ionization 

In this section, we focus on different phenomena which might occur as a result of the 
interaction of the hydrogen molecule with the intense laser fields such as dissociation and 
serious kind of ionization process. Due to the exclusion Pauli principle, the electronic 
wave function must always be antisymmetric. We assume an initial singlet two-electron 
state, so that the spatial part, ip{z\, z 2 , t; Ri(t), R 2 (t)), of the two electron wave function 
is symmetric. As it is schematically shown in Fig. [3j the simulation box is constituted 
by two degrees of freedom which are related to the movement of the electrons along 
polarization directions. One of electrons' direction adjoins with the horizontal axis (z\) 
and the other is on the vertical one (z 2 ). Both of them include the negative and the 
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positive halfs. The symmetry of wave-function in the box is identified by the diagonal 
line (X) that traces from the down left corner to up right corner. This symmetry makes 
some opportunity for saving the necessary CPU time and memory in simulation and 
also increases the rate of the calculations. So, during computation, we can just consider 
the upside part of the simulation box with respect to the X line in the Fig. |3j According 
to indistinguishably, the behaviour of the downside of the box will be same as to the 
upside. For example, ip(zi, z 2 ,t; Ri(t), R 2 (t)) of the H^II) regions in the upside part 
is equal to i/j(z 2 , Zi, t; Ri(t), -^2(0) °f the H^II) for the downside part. Therefore, we 
calculate wave function for the upside part and derive the downside date from upside 
part. 

At the beginning time of the simulation, the hydrogen molecule is placed at an 
initial internuclear distance in the ground electronic state with the nuclei at rest. By 
passing the time, the density distribution of the electron clouds and also the position of 
the nuclei will be affected by the radiation of the laser field which leads to dissociation, 
ionization, and accomplishment phenomena. As long as the system is not ionized and 
dissociated, the population remains in the H2 region of the Fig. |3j 

We now review some different kinds of outgoing of the electrons from region H 2 . 
If due to the irradiation, just one of the two bounded electrons becomes affected; that 
electron goes away from the nuclei and the other one rests under the Coulomb attractive 
forces. As a result of this process, the system goes out from the H2 region and enters the 
neighbouring region. This process is called single ionization. The single ionization 
can be occurred through four pathways which are shown in Fig.[3]by 1-4 arrows. Because 
of the indistinguishability, there are only two distinct pathways for the single ionization 
that can be presented by the arrows 1 and 2 or by the arrows 3 and 4. The pathways that 
are shown by arrows 3 and 4 are completely equivalent to arrows 1 and 2 respectively. 
In the pathway 1, the electron 2 (e 2 ) goes away from nuclei in the positive direction 
while in the pathway 2 the electron 1 goes away in the negative direction which is in the 
opposite direction of the electron 2. The distinguishable movement in the pathway 1 
and 2 is due to the symmetry breaking of the system by irradiation. Therefore, there are 
two distinguishable regions labelled I and II which are related to the region. In the 
absence of the laser field, the initial symmetric wavefunction remains symmetric with 
respect of the negative and positive half. Beyond the H2 region, the force of the external 
field overcomes the Coulomb forces on the ionized electrons. Therefore, outside the H 2 
region, the movement of the ionized electrons is controlled mainly by the laser field. In 
the presence of the laser field after first ionization, it is expected that the next ionization 
take places. In this condition, the distance of the non-ionized electron increases and as 
a result the system goes out of the H^I) and H^" (II) regions which are related to the 
single ionization and enters the neighbouring H 2 + regions. The pathways of the second 
ionization are represented in the upside region with respect to the X line and by arrows 5- 
8 in Fig. [3j Therefore, there are four distinguishable pathways for the second ionization. 
This procedure leads to the sequential double ionization. Figure [3] shows three recognized 
different regions which are related to the double ionization, namely H 2 f+ (I), H2~ + (II), 
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Figure 3. (Color online) Schematic representation of the simulation box used in 
this work. It is constituted by two degrees of freedom pertinent to movement of the 
electrons along polarization direction. One of electrons' direction adjoins with the 
horizontal axis (zi) and the other is on the vertical one (z^). The simulation box are 
divided to some different regions related to the H2, H^, and H^ + . The symmetry of 
wave- function in the box is identified by the diagonal line (X). So, during computation, 
we can just consider the upside part of the simulation box with respect to the X line. 
According to indistinguishably, the behaviour of the downside of the box will be same 
as to the upside. The outer red regions surrounding the inner regions of the simulation 
box show the boundary absorption which lets absorption of the outgoing single and 
double ionization wavef unctions. 

and H2~ + (III). The distinction of the H2~ + (I) and H2~ + (II) regions, like the H^(I) and 
H2~(II) regions in Fig. [3J is due to the radiation of the linearly polarized laser field 
that leads to the symmetry break between the left and right of the wavefunction. The 
behaviour of the system in the third region of the second ionization, i.e. the H2~ + (III) 
region, is completely different. Results in [27] (see Fig. 9) show a higher probability 
of finding the system in H2~ + (III) regions than in regions H^ + (I) and H2~ + (II). This 
is due to the tendency of the system for remaining in a more stable region. When the 
system is in the H^~ + (III) region, the two electrons have gone away from the nuclei in 
the opposite directions and as a result the system has the least feasible inter-electronic 
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repulsion. In the H^" + (I) and H^ + (II) regions, the both electrons are in the same positive 
or negative direction. In this situation, the repulsive force between the two electrons 
makes the system more unstable. Therefore, it is expected that the double ionization 
with opposite directions (H^" + (III)) is more probable than the double ionization with 
the same negative (H^ + (I)) or positive (H2 + (II)) directions. As shown in Fig. [3} the 
boundary absorption has been placed in the end of the regions of the single and double 
ionization which lets absorption of the outgoing wavefunction. 

The indistinguishability of two electrons remains unchanged in the intense laser field 
but the ionization process occurs via various distinguishable pathways. There are two 
different pathways for the single ionization and four different pathways for the double 
ionization. In the double ionization case, two electrons may go out from the right or 
left side of the simulation box. In the third pathway, the two electrons go out in the 
opposite sides. 

IV. Results and Discussion 

The electronic ground state of the hydrogen molecule is calculated by the imaginary 
time propagation method. The obtained potential energy curve is presented in Fig. [TJ 
The electronic ground state energy is zero for the internuclear distance of the 0.45216 
a.u.. The minimum energy is -1.39 a.u. and the related equilibrium distance for the 
classical nuclear motion is 2.13 a.u.. 

At first we study the field free behaviour of the system. The details of the size of the 
simulation box for the field free calculation is as follows. The grid points for each zi and 
Z2 coordinates are 1110. The finest grid size value in an adaptive grid scheme is equal 0.2 
a.u. for both z x and z 2 coordinates. The grids extend up to z x and z 2 =±114 a.u.. The 
size of the absorber regions equals ±14 a.u.. Therefore, the size of the simulation box 
equals 200 x 200, regardless of the absorber regions. In the present laser field, different 
size of the simulation boxes is used for the simulation box that is characterized. 

In the field free simulations, at the initial time, the system is released at the 
different internuclear distances, inner or outer turning points. The calculation shows 
three different sections on the electronic ground state energy curve. The first section is 
related to the internuclear distance smaller than 1.4 a.u. (Ro <1.4 a.u.), the next section 
is related to the interval between 1.4-5 a.u. (1.4 a.u.<Ro <5 a.u.) and the last section 
is for the internuclear distances which are greater than the 5 a.u. (Ro >5 a.u.). Since 
in the first section (R <1.4 a.u.), the nuclei are very close and the Coulomb repulsive 
force is so large, the system is unstable and the dissociation takes place easily. The 
smaller the internuclear distance, the faster dissociation process takes place. However 
the reduction in the population of the electrons and the resultant ionization value is 
ignorable in this section. If the system is released in a internuclear distances belonging 
to the second section (1.4 a.u.<Ro <5 a.u.), which we call the potential well section, 
the ionization and dissociation do not take place which is shown in the Fig. |4j By 
time evolution, the system just oscillates around the equilibrium internuclear distance. 
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More close the initial internuclear distance to the equilibrium distance, the period of 
the oscillation decreases and as a result, the symmetry in the oscillation curve increases. 
At the equilibrium point, there is no oscillation. The antisymmetric which is observed 
in the oscillation cycle in the initial internuclear distances which are far away from the 
equilibrium point, is due to the antisymmetric in the shape of the potential barrier walls. 
For example, the behavior of Ro=1.9 a.u. is relatively symmetric about the equilibrium 
but the behavior of 1.5 a.u. is completely antisymmetric about the equilibrium point 
(2.13 a.u.). Fig. [i] shows that the Ro =1.5 and 3.9 a.u. are the related inner and outer 
turning points. In the third section (R>5 a.u.) the Coulomb force is so weak and the 
bounding between the nuclei is not very strong, so the dissociation does not take place 
unless the nuclei move to a greater internuclear distance at the initial time. 
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Figure 4. (Color online) The time dependent internuclear distances belong to the 
second section (1.4 a.u.<Ro <5 a.u.) in a virtual laser field with a wavelength of 532 
nm and zero intensity. 

In the present of the laser field, at first we discuss the effect of the size of the 
simulation box on the nuclear dynamics. The calculations show that in the higher 
intensity (5 x 10 15 Wcm~ 2 ), the time dependent internuclear distance for different 
size of the simulation boxes is same as 200 x 200. In the lower intensity (5 x 10 14 
Wcm~ 2 ), the results for the various sizes of the simulation box are shown in Fig. [HJ 
This figure shows that during 8 cycles of the laser pulse, there are a small difference in 
the internuclear distances for the various sizes. After 8 cycles, the difference in these 
internuclear distances become considerable and the increase in the size of the simulation 
box does not result in the convergence of the time dependent internuclear distance. 
We changed our explanation in the article as below: We can explain these results as 
follows. In the higher intensity of the laser pulse (5 x 10 15 Wcm~ 2 ), the increasing size 
of the simulation box does not effect on the nuclear dynamics. In this intensity, two 
electrons becomes completely far from nuclei during a few cycles of the laser pulse and 
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a dominate Coulomb explosion occurs between nuclei. Therefore, the increasing size of 
the simulation box does not effect on the behavior of the nuclear dynamics. On the 
other side, at the lower intensity, the electrons do not move completely far from the 
nuclei. Therefore, both the Coulomb repulsion and the population of the electron about 
nuclei determine the time dependent internuclear distance. Therefore, when the size of 
the simulation box increases, the reminded population of the electrons about nucleus 
in simulation box before absorption by boundary absorption becomes slightly more and 
the magnitude of the Coulomb repulsion between nucleus is decreased. 




Figure 5. (Color online) The time dependent internuclear distance for different sizes 
of the simulation box exposed to the laser field with a wavelength of 390 nm and 
5 x 10 14 WcmT 2 intensity. The results for the various sizes of the simulation box form 
200 x 200 to 1000 x 1000 a.u. 2 are shown. 

The behaviour of the system under the action of the laser irradiation is different 
from what is seen in the absence of the laser field. Figure [6] shows the time dependent 
internuclear distance for three different wavelengths and two different intensities. In this 
figure, the size of the simulation box for the 5 x 10 15 Ware 2 intensity is 200 x 200 and 
for the 5 x 10 14 WcmT 2 intensity is 300 x 300. In addition, Fig. [7] represents the time 
dependent behaviour of the population (the residual norm in the simulation box that 
is equal to the total population of different regions in Fig. [3]) in the simulation box for 
three different wavelengths and the two intensities of the 5 x 10 14 WcmT 2 and 5 x 10 15 
WcmT 2 . In this figure, the size of the simulation box for both intensities is 200 x 200. 
The reduction of population is due to the outgoing of the electrons from the absorption 
boundaries, shown in Fig. [3j that results to the single and double ionizations. In a 
separate article, we have shown the details of the single and double ionization [27] . The 
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results of Fig. [6] and Fig. [7] show respectively the magnitude of the internuclear distance 
(dissociation) and ionization for different sections and it is interesting to compare with 
the results of the field free case in Fig. |4| In the first panel in Fig. [6] and [7j i.e. Ro =0.7 
a.u. that belongs to the first section, (R <1.4 a.u.), the dissociation occurs; just like 
what happens in the the absence of the laser field. We also observe a notable decline in 
the norm of the system which means that the electrons' population has gone out of the 
computation box and the ionization has occured. The panels with R = 1.5, 2.13, 3.0, 
and 3.9 a.u. in Fig. [6] and [7] belong to the second section, i.e. (1.4 a.u.<Ro <5 a.u.), 
and in contrast to the field free case both the ionization and dissociation have taken 
place. 

An interesting point in Fig. [6] is that the bond hardening phenomenon does not 
appear anywhere in these classical nuclear dynamics simulations. The bond hardening 
can take places when the laser field causes the nuclei to become closer to each other with 
respect of the field free case. The meaning of the bond hardening in the quantum nuclear 
dynamics investigation that was reported in the previous experimental and theoretical 
researches [HI |9] is vibrational trapping. In the vibrational trapping, the molecular 
wave packet is trapped in a laser-induced potential well. In contrary to the intuitive 
expectations, increasing the laser intensity can lead to the temporary stabilization of 
the molecular bond. However, clear confirmation of bond hardening (or vibrational 
trapping) has remained elusive and might benefit from a fresh look [10J. 

In the last panel with R = 5.0 in Fig. [6] and [7] that is related to the third section 
(Ro >5 a.u.), in contrast to the field free situation, both dissociation and ionization take 
place. In this section, the nuclei go away from each other slower in comparison with the 
second section, but the ionization rate is higher than the second section. 

In all panels of Fig. [6] and [7| we can see that the electrons under the influence of 
the laser field go away from the nuclei which leads to the Coulomb explosion between 
the nuclei. The Coulomb repulsion for the second section is stronger than the third 
section which is due to the smaller initial internuclear distance in the second section 
when Coulomb expulsion is started. 

Figure [6] shows that the speed of the dissociation increase with increasing the 
wavelength or intensity of the laser field. This figure also shows that in the first section 
the wavelength is more effective than the intensity. In fact for the longer wavelengths, the 
nuclei separate from each other more quickly since in this case, the system is irradiated 
for the longer time in each half cycle of the laser field without a change in the direction 
of the laser field. However, in the second section, the influence of the intensity overcomes 
the wavelength influence relatively. In fact in the second section, more the internuclear 
distance is close to the equilibrium point, the intensity increases the speed of the nuclei 
separation more. 

Figure [7] shows that the population reduction starts sooner with increasing the 
wavelength and the intensity. It should be mentioned that due to the increment of the 
stability of the system by getting closer to the initial equilibrium internuclear distance 
(Ro =2.13 a.u.), the population reduction is decreased. This figure also represents that in 
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Figure 6. (Color online) The time dependent intcrnuclear distance for the three 
wavelengths and two different intensities and comparison with the results of the field- 
free in Fig. |4j 

the second section for the smaller wavelengths, the ionization appears before beginning 
of the dissociation which shows that the laser irradiation makes Coulomb explosion by 
forcing the electrons' population to get away from the nuclei. 
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Figure 7. (Color online) The time dependent of the population for the different 
initial internuclear distance for the three wavelengths and two different intensities. 

V. Summary 

In this work and the following article [27], we tried to fundamentally pay attention 
to this subject that how two electrons behave and interact in a two-electron molecule 
with the simplest classical manner of nuclear dynamics. In one part of this article, the 
distinguishable pathways and regions of the single and sequential double ionizations 
were determined and discussed. It is shown that the single ionization can occur through 
four pathways, but there are just two distinct pathways for the single ionization because 
of the indistinguishability of the electrons. This distinguishablity of the two pathways 
is due to the symmetry break between the right and left coordinates which is caused 
by the linearly polarized laser field. In the absence of the laser field, there is not 
such distinguishable pathways for electrons. The result leads to the appearance of two 
distinguishable regions for the region in a linearly polarized intense laser pulse. It is 
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also shown that there are four distinct pathways for the ionization of the second electron. 
We showed that there are two different regions which are related to the single ionization 
and three different ones that are related to the double ionization. Among the regions 
of the second ionization, the H^ + (lll) region is more stable than the two other regions. 
In this region, the two electrons go away from the nuclei in the opposite directions and 
as a result the system has the least feasible inter-electronic repulsion. However in the 
H2" + (I) and H2~ + (II) regions, the two electrons go away in the same positive or negative 
direction. 

In the second part of this article, the time dependent Schrodinger and Newton 
equations are solved simultaneously for the electrons and the nuclei of H 2 respectively. 
These calculations show that there are three different sections on the electronic ground 
state energy curve in the absence of the laser field. In the first section (R <1.4 a.u.), 
the system is unstable and the dissociation takes place easily. However, the ionization 
value is ignorable in this section. In the potential well section (1.4 a.u.<R <5 a.u.) 
the ionization and dissociation do not take place. The system just oscillates around 
the equilibrium internuclear distance. More close the initial internuclear distance to the 
equilibrium distance, the period of the oscillation decreases and as a result, the symmetry 
in the oscillation curve increases. At the equilibrium point, there is no oscillation. In 
the third section (R>5 a.u.), the dissociation does not take place unless the nuclei move 
to a greater internuclear distance at the initial time. In the present of the laser field, in 
all three section, the reduction of the electrons population is appeared in contrast to the 
field free case. The speed of the dissociation increase with increasing the wavelength 
or intensity of the laser field. In the first section the wavelength is more effective 
than the intensity. In the second section, the influence of the intensity overcomes the 
wavelength influence. In the third section (R >5 a.u.), the nuclei go away from each 
other slower in comparison with the second section, but the ionization rate is higher than 
the second section. In the present of the laser field, the time dependent nuclear distance 
is independent of the size of the simulation box at the higher intensity (5 x 10 15 Wcm~ 2 ) 
but at the lower intensity (5 x 10 14 Wcm~ 2 ), the final magnitude of the nuclear distance 
depends on the size of the simulation box. Finally, the bond hardening phenomenon 
does not appear anywhere in these classical nuclear dynamics simulations. 
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